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The classical theory of diffraction at 
plane apertures illuminated by normally 
incident light is modified so that 
diffraction on the source side of the screen 
is taken into consideration and the energy 
transport across the aperture plane is 
described by continuous functions. 
The modified field expressions involve 
the sums and differences of the Rayleigh- 
Sommerfeld diffraction integrals as 
descriptors of a bidirectional flow of 
energy in the near zones on either side of 
the aperture. The theory is valid for 
unpolarized fields, and a pragmatic 
argument is presented that it is appHcable 
to metaUic as well as black screens. The 
modified field expressions are used for 
numerical near-field computations of the 



diffraction profiles and transmission 
coefficients of circular apertures and slits. 
In the mid zone the modified theory is 
reduced to the Fresnel approximation, 
and here the latter may be used with 
confidence. 
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1. Introduction 

This is a continuation of previous papers [1,2] in 
which the physical significance of the classical 
Rayleigh-Sommerfeld and Kirchhoff diffraction 
integrals was assessed and their suitability for compu- 
tations in the near zone was analyzed. The need for 
such computations arises, for example, in the evalua- 
tion of radiometric diffraction errors, where it is neces- 
sary to know the transmission coefficients of the aper- 
tures used for the measurements. The computation of 
these coefficients is a near-zone task even for large 
aperture-detector distances. 

The specific situation considered is a plane aperture 
A contained in an infinitesimally thin screen S that 
occupies the xy-plane of a Cartesian coordinate system 
and is illuminated from the half space z < by a 
normally incident monochromatic plane wave with 



irradiance Eq and wavelength X, The resulting optical 
field is denoted by a scalar wave function, 

(1) 



U{V) = JE,u{n\u{^)\<\. 



and is expressed in the Rayleigh-Sommerfeld theory in 
terms of the surface integrals, 

JkQP 



<(P) = ~JdQ 



In i QP 



, z>0. 



(2a) 






where a metallic screen illuminated hyp- or ^--polarized 
light is assumed.^ The corresponding expression in 



See comments in Sec. 4, below. 
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Kirchhoff s theory, which is usually associated with 
black screens, is 



1 



4k ^ dz QP 

^\ [<'(P) + «Rs(P)L z>o. 



(2c) 



In these equations, A denotes the aperture area, 
P = ^^y-,"^ is the point of observation, Q = (^,?7,0) is a 
point inside the aperture, QP is the distance between 
these points, dQ is the surface element at Q, k=2KlX 
is the circular wave number, and the time dependence 
of the field is assumed as e"'^^ 

Equations (2a,b) were reduced in Ref [1] to previ- 
ously unknown single integrals for the respective cases 
of circular apertures and apertures bounded by straight 
lines, and these were used for numerical computations 
ofw|fs^ andwj^s that involved no simphfying assump- 
tions and could be performed for arbitrarily small 
distances z from these apertures. The numerical 
results obtained were everywhere finite, free of singu- 
larities, and confirmed the well-known prediction that 
u^l^ and w^^s reproduce the boundary values assumed 
in their derivation {du^l^ldz ^\k and u^H -^\ as 
z^O) but not the compatible values (w^g^^ ^ ^^^ 
u Rg/^z -^ ik) which are implied in the classical postu- 
late that the aperture field is the same as the un- 
perturbed geometrical field incident on the screen. 
These inconsistencies obscured the differences between 
the Rayleigh-Sommerfeld and Kirchhoff integrals in 
the immediate proximity of the screen and made it 
impossible to assess their physical significance without 
additional considerations. 

This impasse was overcome in Ref [2] by evaluating 
Eqs. (2a,b) for the special case of a diffracting half 
plane and comparing them to the corresponding values, 
u^ and u^] given by Sommerfeld's rigorous theory 
of half-plane diffraction [3,4]. The agreement was 
remarkably good on the positive side of the screen, 
where the differences (w^g'^^ - Wg^'^^ ) and their derivatives 
were negligibly small even at sub-wavelength distances 
z. Thus, it was decided that the aperture values given by 
the Rayleigh-Sommerfeld integrals are consistent with 
Sommerfeld's rigorous theory, so that attempting to 
improve them would be pointless. 

Accordingly, it became apparent that the real 
problem with the Rayleigh-Sommerfeld and Kirchhoff 
theories was not their failure to reproduce the assumed 
boundary values but these boundary conditions them- 
selves. The classical theories involve "inclination 



factors" which explicitly preclude a backward motion 
of diffracted light and, thus, any perturbation of the 
geometrical field on the source side. On the other hand, 
Sommerfeld's rigorous theory showed that the incident 
light is modified by diffraction before it reaches the 
screen, and therefore the notion of an unperturbed 
incident field is abandoned in this paper by adding a 
diffraction term to the geometrical field on the source 
side. 

The comparison with Sommerfeld's theory also sug- 
gested the need for a further modification of the classi- 
cal theory. The optical field specified by the rigorous 
theory is expressed in the form w^^'^^= w^+w^, where 
u^s''^ obey the same boundary conditions as the 
Rayleigh-Sommerfeld integrals and their components 
Ws and Wg propagate in the opposite directions of the 
incident field and its reflection from the screen and are 
mutually incoherent. In this paper, the Rayleigh- 
Sommerfeld integrals will likewise be resolved into for- 
ward and reverse components defined by 

Wk =-(4s^ +4s) . K =-(4s^ -4s) . ^>0, (3a) 



where the subscript K is used because the forward wave 
function on the left-hand side of this equation happens 
to be the same as the Kirchhoff diffraction integral, 
Eq. (2c). The effective, time-averaged flow of field 
energy is then given by the squared moduli of these 
functions, so that the mutually incoherent, forward and 
reverse irradiances incident on any given area element 
dxdj are given by 



E = E^\u^ 



E = E^ 



^ 



z>0. 



(3b) 



Finally, these quantities will be extended into the 
source space by matching functions so that the overall 
field is continuously differentiable^ in the aperture 
plane and the bidirectional transport of energy through 
the aperture is also expressed by continuous functions. 
The modified theory presented in this paper is valid for 
normally incident light but can easily be adapted for 
obhque angles of incidence. As will be shown, it 
becomes indistinguishable from the usual Fresnel 
approximation in the mid zone z » A, and here the 
latter can be used with confidence. 



That is, continuous witli continuous first derivatives. Tlie purpose 
of tlie derivations in tliis paper is to make the overall field continu- 
ously differentiable with respect to z. Continuity with respect to 
X and y is assured as the diffracted components of the field will be 
expressed in terms of Eqs. (2a-c) and obey the wave equation. 
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2. Modified Field Expressions 
2.1 Derivations 

In addressing the problem of diffraction on the 
source side of a plane metallic screen illuminated 
by normally incident parallel light, it is frequently 
assumed that 



,,(.p>^) — -,,(P'^) 



u\'^'\ z>0, 



Ap,s) 



(4a) 



, ^<0, 



where v^^'^^ is the total field, u^^''^ denotes the field 
components due to diffraction, and e ±e is the 
unperturbed geometrical field on the source side. These 
assumptions appeared first in Rayleigh's papers \_5,6\ 
on diffraction by infinitesimally small apertures and 
show that a continuously differentiable solution for 
yip,s) Yausi obey the boundary conditions 



Xp) ■ 



,(^) 



■l + u 



Xs) 



(p) 



du^^^ du_ 



ip) 



, ^ = 0, 



du 



is) 



dz 



■■2\k + 



du 



(^) 



dz 



(4b) 



(4c) 



These conditions were used by Rayleigh to derive the 
initial terms of Taylor expansions for u ^^''^ for slits 
and circular apertures with dimensions smaller than the 
wavelength of light. Additional higher-order terms 
were calculated by Sommerfeld [4], Bouwkamp [7], 
and others. 

As mentioned above, the Rayleigh-Sommerfeld 
integrals, Eq. (2a,b), will be retained in this paper by 
assuming 



i\'-'=u'l^\x,y,z), z>0, 



(5a) 



and then the second condition in Eq. (4b) and the first 
condition in Eq. (4c) will be satisfied by also assuming 

u['^=-u^ii{x,y,-z\ wi^^=4^(^j,-z), z<0. (5b) 

However the two remaining conditions in Eqs. (4a,b) 
are still not satisfied, so that 3v^V3z and v*^"^ will sfill 
be discontinuous in the aperture plane. 

This failure of Eqs. (4a) can be attributed to the fact 
that the Rayleigh-Sommerfeld integrals are composite 
quantities which can be resolved into the forward and 
reverse field components u^ and u^ in Eqs. (3a); that is, 
^Rs^=^K+^K ^^^ wjfs^= w^-w^. There are no physical 



reasons why these sums and differences should be 
continuously differentiable for z = 0, but on the other 
hand this must be required of Wk and Wk in order to 
correctly account for a continuous transport of energy 
through the aperture. To satisfy this requirement, 
we retain Eq. (5 a) but reverse the signs of e""^ and 
u^^^ in Eqs. (5b), so that 



Xp) 



ti^^(x,j,-z), ,,i^)=-4>(x,3;,-z), z<0. (6a) 



Hence, by applying Eqs. (3a) and letting v = ^(v^^^+v^'^), 



1 {,,{p) 



(v^^^-v^^O, 



v= u^{x,y,z) , z>0 , 

= e^ +u^{x,y,-z) , z<0 , 

v= u^{x,y,z) , z>0 , 

= -e"'^+WK(x,j,-z) , z<0. 



(6b) 
(6c) 



Now, it may be recalled that the Rayleigh-Sommerfeld 
integrals obey the boundary conditions assumed in their 
derivation; that is, 

wj,^^^i)t^^^^ = linA, ^Oon S, z->0, (6d) 
dz 

and hence it follows at once that the scalar field speci- 
fied by Eqs. (6b,c) is continuously differentiable in the 
aperture plane. Likewise, the corresponding forward 
and reverse irradiances, 



E = \u^{x,y,z)f , z>0 , 
^W"^ +u^{x,y,-z)\\ z<0, 

E = \u^{x,y,z)\^, z>0 , 
= |-e"^+w^(x,j,-z)|\ z<0 , 



(7a) 



(7b) 



are continuously differentiable, thus implying a smooth 
bidirectional transport of energy through the aperture. 

Equations (6b,c) and (7a,b) represent the key find- 
ings of this paper. It should be noted that in these 
expressions the roles of u^ and u^ are reversed on 
opposite sides of the screen. That is, u^ appears in the 
expressions for the forward field quantities u^ and E, 
and vice versa. The general properties of these 
modified field quantities can readily be predicted 
from the results reported in Ref [1]; namely that 
the differences between the Rayleigh-Sommerfeld 
integrals u^^^ and u^l are pronounced only in the 
immediate proximity of the screen and vanish in the 
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Fresnel approximation. Thus, a bidirectional exchange 
of energy between the positive and negative sides of the 
aperture occurs only in the near zone, long as WkCx^j^^IzI) 
is appreciably different from zero for values of \z\ on 
the order of a few wavelengths. In the Fresnel limit on 
either side of the aperture plane (|^| » X) the forward 
and reverse fields are unidirectional, and the forward 
field is reduced to the standard expressions in terms of 
Fresnel integrals for slits and Lommel fiinctions for 
circular apertures [8] for z > 0, and to the unperturbed 
geometrical field for z < 0. Similarly, the reverse 
Fresnel field is zero for z > and equal to a Fresnel 
diffraction pattern superimposed on the reflected 
geometrical field for z < 0. 

2.2 Numerical Examples 
2.21 Slits 

As an illustration of the behavior of the modified 
field expression on both sides of the aperture plane. 
Figs, la-c show the forward irradiance profiles 
\_E{x,z) vs x/wf for a slit of width 2w= 10 A and for 
varying distances ±z from the aperture plane. The 
numerical values shown in these figures were comput- 
ed using Eqs. (3a) and (7a) in conjunction with the 
expressions for w^^/^ derived in Sec. 3.3 of Ref [1]. 

Figure la shows that, for z = ± 0.01 A, the modified 
field irradiances are manifestly continuous inside the 
aperture {xlw < 1) and that a modulation of the incident 
field by diffraction also occurs on the opaque portion of 
the screen {xlw < 1). For z = ± A, shown in Fig. lb, the 
diffraction profile on the positive side of the screen is 
already significantly altered in that more light is 
spreading into the shadow, whereas on the negative 
side the modulation of the field is diminished. Finally, 
for z = ± 1 A as shown in Fig. 1 c, the profile on the 
positive side is similar but not yet equal to the Fresnel 
approximation (F, shown as a dashed line) and the mod- 
ulation of the incident field on the negative side is very 
small. This confirms the expectation that the modified 
theory affects only the positive and negative near zones 
in which the Fresnel approximation does not apply. 
For the slit width assumed here, it is estimated that 
the Fresnel limit is reached, within 1 % or better, for 
|z| = 100A. 



2.22 Circular Apertures 

The numerical data presented in Figs. 2 and 3 illus- 
trate the bidirectionality of the field in the positive and 
negative near zones of a circular aperture of width 
2w= 10 A. The data were computed using the mathe- 
matical expressions derived in Sec. 3.2 of Ref [1]. 

Figures 2a,b show the forward and reverse axial irradi- 
ances E{0,z) and E{0,z) for the range -10 A < z < 10 A. 
as given by the closed expressions 



^(0,z)=l-(l+-|) cos[KW-z)-\+-(\+^ 
W 4 



^' 



E{0,z)=-{\-—Y , W = yly^-^ 



W 



(8a) 



(8b) 



which are valid for positive and negative values of z 
and follow readily from Eqs. (9a,b) of Ref [1] and 
Eqs. (7a,b), above. In Fig. 2a, the modification of the 
incident geometrical field is evidenced by the onset 
of pronounced oscillations of the forward irradiance 
E for z < 0. In Fig. 2b, the small but finite values of 
E for z < demonstrate again that the energy flow is 
bidirectional and the reverse field reaches into the 
positive near zone. 

It will be noticed that the reverse axial irradiance 
E{Q,z) in Fig. 2b exhibits no oscillations with respect to 
z. This is due to the fact, illustrated in Figs. 3a and b, 
that E{x,z) always has a maximum for x = 0. On the 
source side, this maximum lies in the reflection shadow 
and is much smaller than the main diffraction pattern 
formed in the region xlw > 1 . 

3, Transmission Coefficients 

The transmission coefficient t of a diffracting aper- 
ture is defined as the radiant flux transmitted into the 
positive half space, divided by the radiant flux incident 
upon it in the limit of geometrical optics. Thus, for a 
plane aperture of area A and normally incident parallel 
light of unit irradiance, 
]_ 
A 



T=^ JdQ ^(Q), 



(9a) 



Because of the two-dimensional nature of tlie diffraction patterns 
discussed in the remainder of this paper, the _y-coordinate will be 
omitted from here on. 



and for two-dimensional apertures of width 2 w which 
are centered on the coordinate origin, as discussed in 
this paper, this is fiirther reduced to 

T = ^}d^ |«k(^,0)|\ (9b) 
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1 ^^^ 1.5 




Q.5 



1 x/w 



1.5 








0,5 



1 x/w 



1.5 



Fig. 1. Forward irradiance profiles on opposite sides of the aperture plane for a 
slit of width 2 w = 10 A. (a) z = ± 0.01 1, (b) z = ± X, (c) z = ± 10 A, F = Fresnel 
approximation. 
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Circuiar Apen^ire 




lAAAA/V 



-10 



z/w 



10 




Fig, 2, Variation of axial irradiance on opposite sides of the aperture for a circular 
aperture of diameter 2 w = 10 A. (a) Forward irradiance E, (b) reverse irradiance £. 



where Eq. (3b) was used and the integral over r\ was 
evaluated as 2 w. According to Eq. (3 a) one finds 



4w i 



lim[<(|,z)]+l 



z>0. 



(9c) 



Here, the aperture value of w^^ was substituted 
from Eq. (6d) in order to avoid computational problems 
that would otherwise arise from singularities for very 
small values of z. The computation of u S^l involves 
lesser singularities and could be performed reliably 
down to z = 0.0003 X. Trial computations indicated that 
the limiting value of T defined by Eq. (9c) was reached 
at the 0. 1 % level for z < 0.003 /L, and consequently the 
results presented below were computed for z = 0.001 X. 



The numerical results thus obtained for the trans- 
mission coefficients of circular apertures and slits are 
shown in Fig. 4 for the range < ^ < 5 71. In both 
cases, these transmission coefficients approach the 
limit T = for ^ = 0, and for larger values of kw they 
exhibit a damped oscillatory behavior. In the case of 
circular apertures the extremes of T occur near kw = Ti, 
2k, ..., whereas for slits they are less pronounced 
and occur near kw = 0.55 71, 1.1 71, ...For large values 
of kw outside the range shown in Fig. 4 both approach 
the limit T= 1, and fiirther computations showed that 
near ^ = 1 00 71 the oscillations of T are still on the 
order of 1 % for circular apertures and less than 0.1 % 
for slits. 
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0.5 



1 x/w 
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Fig. 3. Reverse irradiance profiles on opposite sides of the aperture plane for a 
circular aperture of diameter 2w= 10 A. (a) z = A, (b) z = -A. 
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12 3 kw/n 4 5 

Fig. 4. Transmission coefficients of r circular apertures and slits vs kw. 
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These results can be compared to a large number of 
data that have appeared in the earlier literature for the 
limiting case of very small apertures, kw ^0. In this 
limit the transmission coefficients shown in Fig. 4 for 
circular apertures are superficially similar to, but not 
the same as those computed by Levine and Schwinger 
[9,10]. In the case of narrow slits the results obtained 
here do not at all agree with those published by 
Bouwkamp [7]. These discrepancies will be addressed 
in a subsequent publication. 

4. Concluding Remarks 

There are two aspects of the modified theory present- 
ed in Sec. 2 that deserve further comments: its failure 
to account for polarization effects and the appearance 
of Kirchhoff's integral in the context of a theory in 
which metallic screens are assumed. 

This work was begun in the anticipation that, 
because of their pseudo-vectorial nature, the Rayleigh- 
Sommerfeld integrals could be used to analyze the 
polarization of diffracted light. This anticipation did 
not materialize. For example, the expressions derived 
in Ref. [1] for the axial values of u^^^ and u^^l per- 
taining to a circular aperture illuminated by normally 
incident light differed from each other, although in this 
case symmetry would dictate the absence of polariza- 
tion effects. It also seemed odd that the computed 
values of u^^^ and u^^l were consistently different in 
the near zone but the same in the mid zone, without any 
indication how the degree of polarization could change 
during the free-space propagation of light. Thus it 
appeared that, in spite of the assumption of different 
boundary conditions for u^^^ andw|^\the "polariza- 
tion" effects predicted by the Rayleigh-Sommerfeld 
integrals were implausible. 

Whereas the analysis of polarization effects was 
clearly an objective of Rayleigh's work [5,6], there is 
no indication in Sommerfeld's writings that he had the 
same goal. In his derivation of Eqs. (2a,b) in Ref [4], 
he did not mention polarization at all but stated that 
separate wave functions and boundary conditions were 
required to overcome a well-known mathematical 
inconsistency of Kirchhoff's theory. In his half-plane 
work [3,4], he took the additional step of expressing the 
forward and reverse wave functions in terms of the 
sums and differences of these separate wave functions, 
thus negating any semblance with a vectorial theory as 
these expressions would otherwise imply the inter- 
ference of mutually orthogonal states. Likewise, any 
association of the Rayleigh-Sommerfeld integrals with 



polarized light is negated in the present paper by the 
introduction of Eqs. (3a) in Sec. 1. Accordingly, the 
appearance of Kirchhoff's integral in the modified 
theory has no significance apart from the fact that it 
happens to be the arithmetic mean of u ^^^ and u ^g . 

Given the fact that the modified theory no longer 
pertains to polarized light, a question arises whether it 
is still limited to metallic screens. It may be observed 
that the Eqs. (6b) and (7a) for the forward field are in 
no way altered if the corresponding expressions for the 
reverse field are simply ignored, as if the screen were 
"black." Thus, these expressions might also be useful 
to describe the forward field v produced by a black 
screen, and similarly it might be possible to describe 
the diffraction by partially reflecting screens by simply 
multiplying the reverse field v by a suitable amplitude 
reflectance. These ideas are akin to earlier suggestions 
to define blackness as the absence of reflection [11], 
but cannot be justified theoretically as a metallic 
screen was assumed in the first place. However, from 
a pragmatic point of view it appears that the results 
obtained in this manner would not be far off, and in this 
context we recall Sommerfeld's comment [4] that "a 
slit scratched in a piece of tin foil produces the same 
diffraction pattern, no matter if it is shiny or has been 
blackened^ 
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